function [L, D] = getLandD(altitude, V, alpha, S)
    global mass g
    [rho, sos] = getAtmosphere(altitude);
    Ma = V/sos;
    CL = 0.4172 + 19.41*alpha + 10.17*alpha^2 - Ma*(0.1004 + 0.7536*alpha);
    CD = 0.3042 + 0.02988*CL^2; 
    q = 0.5*rho*V^2;
    L = q*S*CL;
    D = q*S*CD;
end